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ABSTRACT 

Several complex organic molecules are routinely detected in high abundances to- 
wards hot cores and hot corinos. For many of them, their paths of formation in space 
are uncertain, as gas phase reactions alone seem to be insufficient. 

In this paper, we investigate a possible solid-phase route of formation for methyl 
formate (HCOOCH3). We use a chemical model updated with recent results from an 
experiment where simulated grain surfaces were irradiated with 200 keV protons at 
16 K, to simulate the effects of cosmic ray irradiation on grain surfaces. 

We find that this model may be sufficient to reproduce the observed methyl for- 
mate in dark clouds, but not that found in hot cores and corinos. 

Key words: astrochemistry - stars:formation - ISM:abundances - ISM:molecules. 
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1 INTRODUCTION 

Methyl formate (HCOOCH3) is the simplest example of an 
ester. It is derived from the formic acid (HCOOH), where 
a methyl group is attached to the carboxyl group. It is an 
important organic complex molecule that was first detected 
by Brown et al. (1975) towards the Sgr B2(N), the richest 
molecular source in the galaxy located in the Galactic Centre 
giant cloud Sgr B2. It is considered to play a key role in un- 
derstanding the origin of life because it leads to the synthesis 
of bio-polymers. Methyl formate has two structural isomers, 
glycolaldehyde (HCOCH 2 OH) and acetic acid (CH3COOH), 
but it has been reported that this molecule is the most abun- 
dant among these isomers (Hollis et al. 2001). Particularly, 
its column density in the Orion Hot Core was derived to be 
9.4 x 10 15 cm" 2 by Ikeda et al. (2001) and this value was con- 
firmed by Sakai et al. (2007). Methyl formate has also been 
detected in the G31.41+0.31 hot molecular core (HMC); its 
column density is observed to be 3.4 x 10 18 cm -2 around 
this region. Cazaux et al. (2003) observed methyl formate 
in the hot core around the protostellar object IRAS 16293- 
2422. Hence methyl formate seems to be ubiquitous in star 
forming regions. 

Several studies have been carried out on methyl formate 
to understand its formation mechanisms, but it is still de- 
bated whether complex organic molecules form on dust icy 
mantles during the cold phase, on grains during the warm- 
up phase or in the gas phase. Despite a large activation 
energy barrier between protonated methanol and formalde- 
hyde (Horn et al. 2004), Garrod et al. (2006) found that 
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the latter route of formation is viable during the icy man- 
tle sublimation phase. Moreover, very recently, Laas et al. 
(2011) have proposed two Fischer esterification (the acid- 
catalysed reaction of a carboxylic acid with an alcohol to 
give an ester) pathways that occur during the warming-up 
phase, involving protonated formic acid and methanol and 
protonated methanol and neutral formic acid, respectively. 
Both reactions have two channels that correlate to cis- trans- 
protonated methyl formate. They emphasised that methanol 
photodissociation branching ratios and warm-up timescales 
influence the relative ratios between these two geometries. 

The possibility that more complex species are formed at 
low temperatures on surfaces has been investigated before 
(e.g. Charnley 1997, 2001; Herbst & van Dishoeck 2009) and 
in particular methyl formate production on dust surfaces 
was first proposed by Herbst (2005): CO, C and O lead to 
the formation of CH3O and HCO radicals, both known as 
methyl formate precursors. Based on this pathway and us- 
ing energetic electrons at 10 K, laboratory experiments were 
performed by Bennet & Kaiser (2007) in order to produce 
methyl formate and estimate the rate coefficients for this re- 
action. Different experiments involving methyl formate were 
also performed by Gerakines et al. (1996) and Oberg et al. 
(2009), which obtained this molecule after UV photolysis 
of pure methanol and CO:CHsOH ice mixtures. Recently, 
a laboratory study by Modica & Palumbo (2010) has sug- 
gested a new solid state route of formation for this molecule. 
By using infrared spectroscopy in the 4400 - 400 cm" 1 range 
for in situ monitoring the sample during the experiments, 
they simulated a cosmic ion irradiation on a binary mixture 
containing CH3OH and CO ice: methyl formate was there- 
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fore released to the gas-phase after sublimation of these icy 
samples. 

In the present work, we include the Modica & Palumbo 
(2010) experiments in our chemical model UCL_CHEM Viti 
et al. (2004a), by extrapolating from the experiments a new 
rate coefficient for the methyl formate formation on grains. 
Our purpose is to investigate whether surface reactions dur- 
ing the cold phase alone can account for the observed abun- 
dances of methyl formate in massive star-forming regions 
and around low-mass stars. In Section 2 we describe the ex- 
perimental procedure; physical and chemical details of the 
chemical model are summarised in Section 3. In Section 4 we 
qualitatively model four different sources, and present our 
findings. 



2 EXPERIMENTS 

Experiments have been performed in the Laboratory of Ex- 
perimental Astrophysics in Catania (Italy). Solid samples 
were prepared and irradiated in a stainless steel vacuum 
chamber where pressure is kept below 10" 7 mbar. The gas 
mixture to be investigated was injected into the chamber 
through a needle valve where it froze onto the substrate (Si 
or KBr) placed in thermal contact with the tail section of a 
cryostat (10-300 K). After deposition the samples were bom- 
barded by 200 keV H+ ions. Ions are obtained from an ion 
implanter interfaced with the vacuum chamber. The beam 
used produces current densities in the range from 0.1 to a 
few fiA cm" 2 in order to avoid macroscopic heating of the 
target. Infrared transmittance spectra of the samples were 
obtained before and after each step of irradiation by a FTIR 
spectrometer (4400-400 cm _1 =2. 27-25 ^m). More details on 
the experimental procedure can be found in Palumbo et al. 
(2004) and in Modica & Palumbo (2010). 

Pure CH3OH and a CO:CH30H mixture were irradi- 
ated at 16 K with 200 keV H + ions. In both cases, after 
irradiation, the IR spectra show several absorption bands 
which testify for the formation of new molecules not present 
in the original sample. Among these, two bands due to CH4 
appear at 3010 cm -1 and 1304 cm -1 (Palumbo et al. 1999); 
CO2 bands appear at 2344 and 660 cm" 1 ( 12 C0 2 ) and 2278 
cm" 1 ( 13 C0 2 ); H 2 CO is detected at 1720 cm" 1 (Hudson & 
Moore 2000); C2H 4 (OH) 2 (ethylene glycol) is observed at 
1090 cm -1 (Hudson & Moore 2000). Moreover a band is ob- 
served near 1160 cm -1 and is attributed to methyl formate 
(HCOOCH3) and a band near 1067 cm" 1 is assigned to gly- 
colaldehyde (HCOCH 2 OH) (Modica & Palumbo 2010). 

Figure 1 shows the ratio between the column density of 
methyl formate and methanol, N(HCOOCH 3 )/N(CH 3 OH), 
as a function of ion fluence (ions cm -2 ), measured after each 
step of irradiation. The experimental data have been fitted 
with an exponential curve 

N(HCOOCH 3 ) /N(CH 3 OH)=A(l-e" CT *) 

where A is the asymptotic value, a is the cross section 
(cm 2 ) and $ is the ion fluence (ions cm" 2 ). 

The fitting procedure gives A=0.0042 and cr=8.6x 10" 15 

cm 2 . 

To apply the laboratory results to the interstellar 
medium conditions we assume a standard ionisation rate 
(( — 3xl0" 17 s" 1 ) and derive the flux F/sm=0.5 ions cm" 2 
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Figure 1. The ratio between the column density of methyl for- 
mate and methanol, N(HCOOCH 3 )/N(CH 3 OH), as a function of 
ion fluence after irradiation of a CO:CH 3 OH mixture at 16 K. 

s" 1 of cosmic ions in the approximation of effective mo- 
noenergetic 1 MeV protons (see Mennella et al. (2003)). 
Fism must be regarded as an effective quantity. It repre- 
sents the equivalent flux of 1 MeV protons, which gives rise 
to the ionisation rate produced by the cosmic-ray spectrum 
if 1 MeV protons were the only source for ionisation. Fur- 
thermore we assume that the cross section scales with the 
stopping power (S, energy loss per unit path length) of im- 
pinging ions. According to SRIM code (Ziegler et al. 2008), 
in the case of protons impinging on a CO:CH 3 OH mixture 
S(200keV protons)=2.9xS(lMeV protons). 

We defined <jism= Axct/2.9 cm 2 , where the product 
Axtr is the initial slope of eq. (1). The formation rate of 
methyl formate is given by 

R(s" 1 )= (j/ SM (cm 2 ) x F/ SM (cm" 2 s" 1 ) 
Using a ism and Fism as estimated above we obtain 
R=6.2xl0" 18 s" 1 . 



3 CHEMICAL MODEL 

The chemical model is UCL.CHEM, developed by Viti & 
Williams (1999) and Viti et al. (2004a). UCL.CHEM is a 
time and depth dependent gas-grain chemical model which 
estimates the fractional abundances (with respect to the to- 
tal number of hydrogen nuclei) of gas and surface species 
in every environment where molecules are present. For this 
purpose it was adapted to model low and high-mass star 
forming regions. The model performs a two-step calculation: 
Phase I starts from a fairly diffuse medium in atomic form 
and undergoes a free-fall collapse until densities typical of 
the gas that will form hot cores and hot corinos are reached 
(10 7 cm" 3 and 10 8 cm" 3 respectively). During this time 
atoms and molecules are depleted onto grain surfaces and 
they hydrogenate when possible. The depletion efficiency is 
determined by the fraction of the gas-phase material that 
is frozen onto the grains. This approach allows a derivation 
of the ice composition by a time-dependent computation of 
the chemical evolution of the gas-dust interaction process. 
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The initial elemental abundances of the main species (such 
as H, He, C, 0, N, Si and Mg) are the main input for the 
chemistry. In our model surface reactions only occur during 
this phase. Phase II is the warming-up phase and follows 
the chemical evolution of the remnant core when the hot 
core itself is formed. It simulates the effect of the presence 
of an infrared source in the centre of the core or in its vicin- 
ity by subjecting the core to an increase in the gas and 
dust temperature. We derived the temperature of the gas 
as a function of the luminosity (and therefore the age) of 
the protostar. The treatment of the evaporation of the code 
is either time-dependent in which mantle species desorb in 
various temperature bands according to the experimental 
results by Collings et al. (2004) or instantaneous in that all 
species will desorb off grain surfaces in the first time-step. 
The temperature profile is another input parameter. For this 
study, the model does not take into consideration the des- 
orption of molecules due to UV photons and cosmic rays. 

Most of the gas-phase reactions that occur in the ISM 
are collected from the UMIST database (Woodall et al. 
2007). Because of the failure of gas-phase chemistry to re- 
produce abundances of particular molecules, in the present 
work a gas-grain network reaction is also taken into account 
and particularly, we focused on methyl formate reactions. 
Table 1 shows on the top the updated gas-phase reaction 
rate coefficients as listed in the most recent database KIDA 
(Wakelam 2009). At the bottom, we report solid state paths 
of formation of methyl formate, taking into account the new 
experimental rate coefficient for this reaction. Methyl for- 
mate is both one of the reactants (R) and product (P); a, 
/?, and 7 are the rate coefficients. The first six are bimolec- 
ular reactions, where an exothermic proton transfer occurs 
between an ion and a neutral species. The last one is an 
electronic recombination, in which the recombination of a 
positive ion and an electron results in the dissociation of 
the molecule. In both cases, a Koiij formula (a modified ver- 
sion of the Arrhenius equation) is adopted to evaluate the 
reaction rate: 

where T is the gas temperature, a is the cosmic ion rate, 
/3 contains the dipolc effect and 7 represents the probabil- 
ity per cosmic ray ionisation, that for the updated methyl 
formate reactions is always zero. /3 includes the enhance- 
ment of ion-neutral rate coefficients for cases in which the 
neutral has a large, permanent electric dipole moment. This 
effect results in a rate coefficient that has a Ti dependence 
at low temperature (10 - 50 K) as observed by Herbst &: 
Leung (1986). In the bottom panel of Table 1 the surface 
path of formation of methyl formate is summarised: one of 
the most abundant species in the ISM is carbon monoxide, 
that undergoes successive hydrogenations in the presence of 
a large amount of hydrogen (Fuchs et al. 2009; Watanabe 
et al. 2007). Consequently, methanol can be formed. After 
cosmic ion irradiation, CO reacts with CH3OH to produce 
methyl formate (see Modica & Palumbo (2010)). The m be- 
fore the molecular formulae has been inserted to underline 
the role of the surface chemistry that acts as a catalyst for 
the reactions with large energy barriers. The new experi- 
mental rate coefficient for this reaction is 6.2 x 10~ 18 s _1 , 
as shown. 



4 RESULTS 

Figure 2 shows the fractional abundances of selected species 
(from Table 1) as a function of time during Phase I. The solid 
lines represent the abundances of methanol, carbon monox- 
ide and methyl formate when the new solid state route is 
included in the model. The dashed line is methyl formate 
obtained from the model computed with the same chem- 
istry as in the original code (Viti et al. 2004a). In the latter 
case, the reactant trends are the same as those shown in the 
updated version of the code. While the reactants, CO and 
CH3OH, reach a fractional abundance of 4 x 10 and 2 x 
10~ 6 respectively, methyl formate is produced in a negligi- 
ble quantity (< 10~ 21 ) onto the grains. By considering the 
solid state path of methyl formate formation with the ex- 
perimental rate coefficient of 6.2 x 10~ 1S s _1 (as described 
in Sections 2 and 3) a value of 1.8 x 10~ 9 is reached for its 
fractional abundance (solid line). In order to compare our 
theoretical abundances with those from the observations, 
we also run Phase II for a typical hot core. In order to fa- 
cilitate the calculation we set an instantaneous evaporation 
instead of a time-dependent sublimation for molecules that 
sublimate from the grain into the gas phase. As a result we 
obtained a detectable abundance of methyl formate with a 
value of 1.7 x 10~ 9 . An alternative route for the synthe- 
sis of methyl formate in a hot core was previously investi- 
gated by Garrod & Herbst (2006), that proposed both gas- 
phase and grain-surface processes to produce most species 
during the warm-up phase. They estimated that the gas- 
phase/accretion path, involving protonated methanol and 
formaldehyde, is responsible for the formation of 25% of the 
total HCOOCH3 present on the grain surface before subli- 
mation. Indeed, during the warm-up phase, the surface re- 
action between CH3O and HCO radicals is allowed and it 
is thought to be the dominant mechanism; when significant 
amount of H2CO begins to evaporate, the methyl formate 
formation on grains is inhibited and formaldehyde is pro- 
tonated to produce methyl formate on gas-phase. These re- 
actions are more efficient in hot corinos than in hot cores. 
They found a trend where methyl formate reaches a frac- 
tional abundance of about 10~ 8 at 10 6 yr. Both mechanisms 
seem therefore to be viable. 

In order to compare the theoretical abundance of 
methyl formate with the observational data, we model two 
high massive forming regions, the Orion Hot Core and the 
G31.41+0.31, a low-mass source, NGC1333- IRAS2 and the 
cold gas near the region Bl-b in Perseus using our two- 
phases UCL_CHEM. The physical and chemical parameters 
for each model are listed below. In Phase II of all models we 
now assume a time-step sublimation instead of an instan- 
taneous evaporation for molecules that sublimate from the 
grains into the gas-phase during the warm-up stage of star 
formation process, as described by Viti et al. (2004a). 

Orion KL Hot Core 

The Orion Molecular Cloud 1 (OMC1) is a dense clump 
at 480 pc that contains several distinct infrared-emitting 
regions. Among these the Kleinmann-Low (KL) region is 
composed of four different components: the hot core, the 
compact ridge, the plateau and the photodissociation re- 
gion (PDR) surrounding the quiescent gas. While the hot 
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Table 1. Gas- phase and solid state routes of methyl formate destruction and formation. The 
top panel shows the gas-phase reactions involving methyl formate, whose rate coefficients are 
taken from the KIDA database. The bottom panel contains the new experimental path of methyl 
formate formation on grain surface as investigated by Modica & Palumbo (2010). 
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Table 2. Observed and predicted column densities towards Orion KL Hot Core, G31.41+0.31, 
NGC1333-IRAS2 and Bl-b core. 



Source 


Type 


Distancc(pc) 


T rot (K) 


Observations (cm 2 ) 


Model(cm -2 ) 


Orion KL Hot Core 


High mass 


480 


250 


a 6.9xl0 16 


3.5xl0 14 


G31.41+0.31 HMC 


High mass 


7900 


300 


6 6.8xl0 18 


9.3xl0 15 


NCG1333-IRAS2 


Low mass 


220 


38 


c 5.8xl0 14 


6.4xl0 12 


Bl-b 


Dark core 


350 


<30 


d 8.3xl0 12 


1.2xl0 12 



a Remijan et al. (2003) 

^ Note that the total methyl formate abundance 
given in Cesaroni et al. (1994) 

c Note that the total methyl formate abundance 
given in Bottinelli et al. (2007) 
d Oberg et al. (2010) 
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Figure 2. Fractional abundances (with respect to the total num- 
ber of hydrogen nuclei) of selected species as a function of time. 



core and the plateau are characterised by elevated temper- 
atures, the ridge consists of an extended and cooler region 
with quiescent material. We focused on the Orion hot core as 
it is one of the richest known astronomical sources of molec- 
ular lines. The main heating mechanism in this location is 
radiative, and gas kinetic temperatures are close to 200 K. 
Using UCL.CHEM Lerate et al. (2008) computed different 
chemical models to simulate each component of the KL re- 



gion. In the present work we used the physical and chemical 
parameters for their best-fitting model. 

A comparison between our model calculations and the 
observational column densities through the Orion Hot Core 
is shown in Table 2. Observational data are taken from 
Remijan et al. (2003), that have surveyed this source with 
the BIMA Array with angular resolution of about 2"-5". 



G31.41+0.31 

G31.41+0.31 is another example of a multiple massive star 
formation region located at 7.9 kpc. It is associated with 
an hot molecular core (HMC) of 0.048 pc in size, where the 
temperature is found to be about 300 K and the luminosity 
is about 3 x 10 5 solar luminosities (Churchwell et al. 1990). 
It is in fact the region where the simplest sugar, glycolalde- 
hyde (and also one of the methyl formate isomers) was first 
detected by Beltran et al. (2009) through IRAM PdBI ob- 
servations. Beltran et al. (2009 also used UCL.CHEM to try 
to reproduce the abundance for glycolaldehyde. We adopted 
the physical and chemical input parameters from their study 
for our model. 

We compared our theoretical results with those given 
by Cesaroni et al. (1994), that performed high resolution 
observations by using IRAM Plateau de Bure Interferome- 
ter. 
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NGC 1333-IRAS2 

Located in the cloud NGC1333, belonging to the Perseus 
complex at 220 pc in distance, IRAS2 is a low-mass binary 
protostellar system, including IRAS2A and IRAS2B, sep- 
arated by 30". An infalling envelope, a circumstellar disk 
and multiple outflows seem to be associated with them. The 
advantage of observing such nearby objects is that these 
sources can be spatially resolved with millimetre interferom- 
eters. As a result we adopted a multipoint model in order 
to reproduce the density profile of the hot corino associated 
with IRAS2, as reported in Maret et al. (2004). We also used 
the temperature profile for hot corino given in Awad et al. 
(20f0). 

We evaluated the total column density for the methyl 
formate as described in Viti et al. (2004b): 

N{HCOOCH 3 ) = ^2(X xA v xL)x N(H 2 ), (2) 

where X is the methyl formate fractional abundances 
and A v is the visual extinction at length L (pc). 

Bl-b core 

The protostar associated with the Bl-b core is believed to 
be between a pre-stellar and Class protostar evolutionary 
stage. Two sources are identified with the Bl-b core, de- 
signed as Bl-bN and Bl-bS, separated by 20" in the north- 
south direction at about 350 pc in distance (Hirano 1999). 
The physical parameters of two sources are very similar 
{Tdust 18 K, mass of 1.6-1.8 solar masses and luminosity 
about 2.6-3.1 solar luminosities). We qualitatively model the 
Bl-b core running our two- phases model; for this core, dur- 
ing Phase II the temperature only rises up to 30 K; this 
implies that only the weakly bound species partly sublimate 
due to thermal desorption (see Collings et al. (2004)); hence 
any methyl formate formed on the grains would remain in 
solid phase. However, non thermal desorption processes are 
known to be efficient even at 10 K (see Roberts et al. (2007); 
Boland & de Jong (1982); Hasegawa & Herbst (1993) and 
Bringa & Johnson (2004)). We have therefore ran Phase I 
including some non thermal desorption mechanisms (as in 
Roberts et al. (2007)). We compare our results for methyl 
formate with those reported in Oberg et al. (2010). Note 
that our computed methyl formate is a lower limit as our 
non thermal desorption mechanisms did not include direct 
UV photodesorption. 

Table 2 reports all our model results and compares the 
theoretical predicted column densities for methyl formate 
with those observed in the sources described above. A rela- 
tively good agreement within a factor of 10 was found for the 
dark core associated with Bl, but the theoretical values are 
too low for the cases of warmer star forming regions. In fact 
around protostars the flux of ions due to stellar flares should 
also be considered. As discussed by Garozzo et al. (2011), the 
effects induced by cosmic rays on icy grain mantles during 
the collapse phase era are comparable to the effects induced 
by stellar flare ions during the warm-up phase. 

In summary, using the rate coefficient derived from lab- 
oratory experiments of ice irradiation, we investigated the 
viability of a cold solid state path to form methyl formate 
during the evolution of protostellar cores when the tempera- 
tures are low (10 K). We find that we cannot reproduce the 



abundance of methyl formate found in hot cores and hot 
corinos without invoking gas- or solid-phase reactions which 
necessitate high temperatures (and that can therefore only 
occur during the warm up phase once the star is formed), 
while we were able to reproduce the observed abundances 
in dark clouds. This work supports the idea that cosmic ion 
irradiation of icy grain mantles may be able to contribute 
to the production of the methyl formate observed in dense 
molecular clouds, but that other routes of formation are also 
required. 
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